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SUMMARY 


A general analysis capable of predicting perfonaance 
characteristics of vertical axis windmills (VAW) has been 
developed. This analysis includes the effects of airfoil 
geometry, su|^ort struts, blade aspect ratio, windmill solidity, 
blade interference and curved flow. The results from the 
analysis are coqsared with available wind tunnel results for 
a catenary blade shape. A theoretical performance curve for 
an aerodynamically efficient straight b'lade conflguraticm is 
also presented. In addition, a linearized analytical solution 
applicable for straight blade configurations is developed. 

This analysis is useful for parameter studies when the 
aerodynamic characteristics of the airfoil can be assiimed 
to behave linearly with angle of attack. The errors introduced 
by this assumption are demonstrated by ccunparlng numerical 
with analytical results. A listing of the computer program 
developed for numerical solutions of the general performance 
equations is included in the appendix. 
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SYMBOLS 



aspect ratio 

swept or projected area in vertical plane 
lift curve slope for two diiaenslonal airfoil 
drag coefficient 
lift coefficient 
Bonent coefficient 

coefficient of the average aonent produced during 
one revolution of the windsill 

e 

power coefficient 

pressure coefficient 

pressure coefficient - lower surface 

pressure coefficient - upper surface 

thrust coefficient 
chord length 
drag 

chordwise distance from airfoil leading edge to 
center of pressure 

height of windmill 

lift curve slope including aspect ratio and downwash 
corrections 

lift 

length 


M 


moment 



SYMBOLS CONTIHUED 


nomal force 
ntober of blades 
power 

total wind stream power 
dynamic pressure 
tip speed to wind speed ratio 
Reynolds nuober 

distance fron axis of rotation to half diord point of 
blade element 

distance from axis of rotation to any particular point 
along blade chord line 

spanwlse distance 

thrust 

wind speed at win^taiill 

wind speed relative to blades 

wind speed in free stream 

chordwlse distance from leading edge 

vertical distance 

angle of actack 

local blade slope angle 

angle in plane of rotation between r and r* directl<ms 
angle of rotation from wind direction 
chordwlse distance from half chord point forward 
air density 
rotational speed 
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STI0OI.S OONTmOED 

Subscripts : 


a 

average 

max 

MUClUUB 

T 

relative 

T 

total 

0 

referent 

OD 

free streea 


Superscripts: 

— nondiaensltmal 


tine derivative 
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IHTRODOCTIOK 

As a ccmsequence of the current shortages of fossil fuels 
interest is being revived In alternate sources of energy. One 
means for converting solar energy to a more useful form is 
through the use of windmills. Windmills may take many forms 
from drag devices to those which utilize the efficiency of 
airfoil shapes. An Interesting survey of various types of 
windmills Is given In reference 1. The most common windmill 
using airfoil shapes is the conventional multib laded horizontal 
axis type. These range from two blades to as high as twenty-seven 
blades depending on application, and the blades may range from 
propeller type to flexible sails. Analyses of horizontal axis 
windmills are well documented. 

Another category of windmill is the vertical axis (VAW) or 
Darrleus windmill. The primary attraction of a VAW is the 
simplicity of manufacture compared to a HAW, and the -manner in 
<Aich the loeds developed during operation are resisted by the 
structure. These differences Indicate that if the efficiency 
of a VAW is equal to or not much less than a conventional 
windmill, use of a VAW might result in a more economical 
energy conversion system. 

In recent years a small body of data has been ptibllshed on 
the vertical axis windodll. Some of these reports present the 
results from tests of the vertical axis concept, ref. 2 and 3, 
while others present various analyses, structural, aerod 3 mamlc. 



and econoolc, of vertical axla windallls, references 4 throt^h 7. 

This report presents a general analysis of vertical axis 
wlndmlUa. Airfoil theory Is used to determine the torque ' 
generated by the blades ^out the axis of rotation, and 
anatentun theory Is used to determine the effective wind velocity 
at the blades. These yield a set of Integral equations which 
require an Iterative approach for their solution. The nonlinear 
characteristics of the blad'^ airfoils are accounted for In the 
general anal 3 rsls. A more simplified analysis assuming linear 

a 

aerod 3 rnamlcs and appUcable to stral^t blade configuration 
Is also presented. 

The results of this analysis In addition to being In gtid 
agreement with available test data on one VAW concept Indicate 
that the VAW can be designed such chat their efficiency Is com- 
parable to the best horizontal axis designs. 

This report presents the details of the general nonlinear 
analysis and the lineariz*id analytical solutions for the straight 
blade concept. The result produced from the general analysis 
are compared with existing test lata. The results from a 
parsa^ter study of VAW's using trie general analysis and the 
linearized solutions are compared and limitation of the linear 
analysis are identified. One straight blade configuration having 
significantly improved aerodynamic performance while maintaining 
the sane structural capability of the tested configuration is also 
analyzed and the results are presented. Also Included in the Appendix of 
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this report it a lltting, of a eonpttttr pri^at for gtatral mtlyata 
of vortical axla vindailla along with input and roaulta for a aanple 
eaaa. 
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GENEBAL TH^IH 

In this section, gecerel nondlmensional equations for 
aerodynaaLc torque and average vlndmlll thmst are developed by 
use of actuator disk theory and blade element aerodynamics. Real 
effects sudt as aspect ratio, dovnvash and flew curvature are 
accounted for. An average nondimensional power coefficient is 
defined and Its relation to wind speed and average windmill 
torque is given. 

Defining Equations 

The basic elements of the vertical axis windmill are 
shown on figure 1. It is con^osed of an arbitrary number of 
blades (the figures are shown with two) attached to a rotating 
shaft which is nominally perpendicular to free stream velocity 
vector. The blades shown in the figure are straight; however, 
this is not a requirement. The blade cross section must be 
such that when subjected to an angle of attack a lift force 
is generated. The performance of the windmill is directly 
related to the ratio of this lift force to the blade drag. 

Since the angle of attack which the blades experience is cyclic 
the blade cross sections are restricted to sjnmnetrlc shapes. 

The coordinate system used in this analysis is shown in figure 2. 

The origin is located at the midpoint of the rotor with the y axis 
corresponding to the spin axis. The wind velocity is arbitrarily 
assumed to be directed along the x-axls. The rotor position as given 
by the angle G is measured from the x-axls. The local radius is measured 
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from tha y-axis and the local slope » 6 » is measured «lth respect to the 
radius vector. 

The determination of local lift and drag characteristics of a blade 
element represents au unsteady flow problem. However, it has been assumed 
that the static aerodynamic data for the subject airfoil section is 
applicable . 

Consider an arbitrary blade section as shown in figure 3. For a 
blade eleaent the Incremental forces acting in radial and tangential 
dlrecticms are derived as follows. The total velocity at any location 

s 

on a blade is, 

- (ru - V sin 0)^ + (V cos 0)^ (1) 


The angle of attack is given by 


tan a ® 


V cos 0 sin B 
rw - V sin 0 


( 2 ) 


The incremental forces developed as a consequence of the angle of 
attack are defined as, dL, in the direction normal to the relative 
velocity vector, and dD, in the direction of the relative velocity. 

These fCirr-s act to produce a moment about the axis of rotation given by 
dll =» (dL sin a - dD cos a) r + (dL cos a + dD sin a) sin 8(-|- - d) (3) 

The elemental lift and drag force can be written as 

dL ■ Cj^ q cds (4a) 

and 


dD • Cp q cds 


(4b) 
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«h«re the eleaent lei^th» de, le given by ds * dy/sln 6. 

Std)8tltuting eq. (4) into eq. (3) yields 

dM ■ (Cj^ q c(dy/sln 6)sin O - Cjj q c(dy/sin $) cos o)r 

+ (Cj^ q c(dy/sin $)cos o + q c(dy/sin B) sin a) sin 6 (y - d) , (5) 

where q - J P« vj • 

Integrating over y yields the total monk>at due to the blade as a 

function of the angle 0. 
h/2 

M(6) - / (C^ .in a - Cj cos o.) 

-h/2 

h/2 

+ / (Cj^ cos Ct + Cjj sin a) qc (j - d) dy (6) 

-h/2 

The following nondiTOnsionalization scheme is used; all velocities 

are nondimenslonalized using the free stream velocity* V^, and all 

lengths are nondimensionallzed using the maximum radius of the 

tfindmlll blades r . The nondimens ional variables are denoted 
max 

b> a ber, i.e, V *» With this notation equation (6) becomes* 


M(0) 




h/2 

-h/2 


- - Tr 2 

r c dy 
sin 8 


P V 

at 


(r )■ 

max' 


h/2 
/ <Ct 


cos ot + Cjj sin ct) c(y - d) dy 


-h/2 


(7) 



u 


If a ■o a cB t coafflciant, (^, ia defined as 


Si 


M 


2 



max 



( 8 ) 


idiere 4 la the area swept out by the wlndnill then eq. (7) becoaes 
8 

h/2 

CgW) . i / «C^ .In a - Cp CO. 

»-h/2 


+ (Cj^ a>s a + Cjj sin a) (-j - d) } c dy 

This can be rewritten as 

h/2 e ^ V 2 

*. 0 

+ - •7) sin 0 (Cj^ cos <* •»■ sin a)] dy (9) 

r c 

Equation (9) holds for one blade, if there are n blades in the 
configuration the contribution for each blade can be obtained 
bj' replacing 0 in equations (1) and (2) by 0 + (i - 1) A0 and 
r( -evaluating equation (9). 

Doing this the total soment coefficient can be written 
i =» n 

\ ^ Si 

^ i - 1 

tdiere A0 » ^ , 

n 
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The integral of the toted o(»ient throu^ one revolatlcm of the wiml^ll 
l8 given by 



n 


2ir/n 




(0) d0 


(U) 


The average iK)rent Is then given by 


\ 

<=M -sr 

a 


2if/n 
f 

2it •' 
o 


(0) d0 

The average potrer developed by the windniU is given by 


( 12 ) 


P = 


aai 


2tr 


2ir/n 

/ Hj, (0) d0 


(13) 


Defining a pover coefficient C as 

P 




vjhere 


P = -I E v j 
t I ^ 


^ PQO A 


represents the total kinetic energy available per unit time in the 
stream tube. Using the definition of eq. (3), eq. (13) can be 
written In terms of the average moment coefficient as 
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P 



thereby gl'viag 


(^) Si 

V_ a 


(14) 


or 


C » R 

P \ 


(15) 


where 


R 




In nondlaensional form equations (1) and (2) become 
V ^ a (R r - V sin 0)^ + cos^ 0 (14) 


and 


V cos 0 sin B 


tan a 


r R - V sin © 


( 17 ) 
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Moaentun Loss «od filndBlll Thrust 

The thrust sctlng on any device extracting energy from the Wind 
Is directly related to tite diange In soeentiai which occurs In the 
wind strean. Using ■oaentuB theory, reference 8 gives 

« 2 p„ V (1 - V) vj (18) 


or 


Cj = 4 V a - V) 

lAere is an average thrust coefficient and is defiled as 
a 



i p. ''J 


A 

s 


(19) 


In terns of the elemental forces acting on a blade section the average 
thrust coefficient can be written 

2ir/n _ h/2 - 2 

S = ^ ^ « + Cy sin a) cos 0 sin 6 

® o -h/2 


+ (Cj^ sin a - cos a) sin 0 ]dy d0 (20) 

The determinatioa of and V for specified values of requires 


an iterative procedure 
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Real Effects 


Although deteniolng the flow field around a rotating windadll airfoil 
constitutes a 3 dlnensional unsteady flow probleia, the asstapticm has 
been aade in this analysis that the flow is steady and that a two 
dlaensional analysis can be used. However, a number of other effects 
including aspect ratio, downwash and curved flow can influence the 
rotor performance and some attempt has been made to evaluate these. 

Aspect Ratio . - During a rotation through 360 degrees a blade 

section will experience large angles of< attack as the ratio of 

ro>/V becomes small. When this occurs blade stall will result and 
00 

to account for stall effects, it is necessary to use experimental 
data on section lift, drag and moment coefficients. These data 
are generally presented for infinite aspect ratio wings with 
appropriate aspect ratio correction being required. From reference 
9 using Prandtl's lifting line theory, aspect ratio correction 
can be written as 


C 


L 



/ (1 + 


a 

o . 
TT AR' 


(21a) 


and 2 

S S n AR 
o 

where ^ represents the lift and drag coefficient for an 


(21b) 


o o 

infinite aspect ratio wing, and AR « h/c. 

Downwash . - For a mu Itib laded configuration each blade is in the 
downwash from the preceding blade thus the angle of attack is reduced by 
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tills dcwivash angle. Again using the spproadi outlined in reference 9 
tile effect of dovnuash can be accounted for by using an aapect ratio 
correction. This yields 

1 - a / ir AR 

^ 1 + a / ir AR ^ 
o o 

Plow Curvature . - Since the airfoil is traveling in a circular arc 
the angle of attadc varies along the length of the chord, nie effect of 
fleer curvature is a functlcm of the spin rate and the chord to radius 
ratio. This effect can be evaluated by considering the flew over 
a flat plate. The force distribution on a plate of unit width in a 
unifoxa stream of velocity is given by 


dN 

dz 




For a flat plate the velocity ratio at angle of attack is given 
by reference 10 as 


(23) 


= cos a ± sin a —) , (24) 

''oo * 

where the (+) and (-) signs refer to the upper and lower surfaces, 
respectively. Now the relation between V/V and C is 

^ ' OO p 

C - 1 - (V/VJ^ (25) 

P 

Substituting yields 

— ■ 4 q cos a sin a —) (26) 

dx ^ X 


Since the blade is actually moving along a curved path, the local 
angle of attadc will vary. Referring to figure 4 the local angle of attack 
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can ba vrittan aa 


cot a 


r’w COB e - sin 6 
00 

V- coe 0 - r*w ain e 


cot a 


rw - Bln 0 

w 

cos 0 - 5 w 

00 


cot a 


R - sin 0 
cos 0 - 1 R 


lAere R * ru)/V_ . 


In term of the coordinate x oeasured from the leading edge of 


the plate this becoaes 


cot a * 


R - sin 0 

cos 0 + R(x - c/2) 


In terns of the lift distribution on the plate eq. (26) becomes 


dL , 2 ^ . ,c - Xx 

B 4 q cos a sin a (—3 — ) . 

dx X 


In the limiting case c o and the expression for angle of attack becomes 


cot o_ 


R - Bln 0 
cos 0 


Defining a reference lift distribution as 


(^) » 4 q cos^ a sin a 

dx X 



Ae effect of curvature can be diaplayed aa a function of the ratio of 
lift dlatrihutlona » «lth and without curvature effecta. That ia 
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(dL/dS) 

(dL/dx) 


2 

cos g aln g 
2 

cos a sin a 
o o 


(30) 


The ratio of lift with and without curvature effects can be written 


/ cos^ a sin a d X 
o 


^2 

f cos^ a sin a d x 
' o o 

o 


/ cos^ a sin a d x 
o 


cos a sin a c 
o o 


(31) 


The nuBwrator can be integrated givit^ 


R - sin 0 
R 


(cos - cos 02 ) 


cos a sin a c 
o o 


(32) 


where 


a 


tan 


-1 


tan 




tan 


( 


cos 6 


-) 


-1 


( 


R - sin 0 

cos 0 - R c 

2 _ 

R - sin 0 

R c 


-1 


cos 0 + 


(- 


R - sin 0 


This correction term has been applied to account for flow curvature 
effecta. For most typical configurations the correction is less than 
one percent. However, if c becomes large this effect would become 
more important. 
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VERIFICATKW OF GEHEKAL THEOBY 

In this sectlim, the difficulties of modeling windmill . 
performance using airfoil data with corrections for three 
dimensional flow effects are discussed. Theoretical results 
are compared to test data. 

Reynolds Number, Roughness and Crossflo'f 

As Indicated In reference 11 the NACA 0012 airfoil falls Into 
the category which experiences either leading edge or trailing edge 
stall depending on Reywilds Nunber. 

Furthermore, the sectlcm characteristics are very dependent on 
surface roughness. When these effects are cosblned with the three 
dimensional nature of the flow over the windmill blades It becomes 
apparent that using two dimensional data at a fixed Reynolds Nunber 
represents a significant approximation. 

Fortunately, use of two dimensional data yields a good 
approximation to the actual performance of a windmill; however, 
differences between measured and computed performance can easily 
be explained within the framework of these aerodynamics uncertainties. 
The data used in this analysis come from reference 12, however, 
some cases are presented in which these data have been modified to 
simulate leading edge stall and surface roughness. 

Comparison with Test Data 

The results predicted by the analysis for a two blade con- 
figuration (n 2) were compared with test data from references 
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3 and 4. The blade ahape la a catenary with maxlmua radlua of 2.13 m 
(7 ft) and a hel^t of 4.27 m (14 ft). The shape of the upper half 
of a blade la represented by the curve shoira on figure S. The test 
Reynolds number was about 0.3 x 10^ consequently low Reynolds nuod>er 
data were used for the NACA 0012 airfoil. The airfoil data were 
obtained from reference 12 » and a tabulation of lift and drag 
coefficient and center of pressure versus angle of attack is presented 
In table 1. These data are for a Reynolds Ifumber of 0.5 x 10^. 

From the theoretical curve It Is seen that the catenary V4W 

a 

does not begin to produce sufficient power to overcome the drag 
due to stall imtll a velocity ratio of about 2.5 is reached. At 
this point the power coefficient increases with spin rate until 
a peak value of about .37 is reached. The power coefficient then 
decreases until at a velocity ratio of R ~ 10 It reaches zero. 

This determines the no load spin rate at which the windmill would operate. 

The test data follows the same trends; however, there are 
differences between these results. The test data seeus to indicate 
that stall occurred at a higher velocity ratio and more abruptly 
than predicted, possibly indicating leading edge stall. Also, 
the no load velocity ratio reached in the tests was slightly greater 
than 6 indicating the zero lift drag was greater than predicted by 
smooth airfoil data. To demonstrate the effect of these variations 
on the predicted performance curve a computation was made using a 
zero lift drag coefficient representative of a rough surface. 
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Increoslog reduces the peak efflclencgr aoBCidiat and causes the 
o 

no load tip speed ratio to drop to a value of about 8.5. To 
doaonstrate the effect of leading edge stall at low Reynolds nind>er 
the lift and drag coefficient data were loodlfied to reflect a 
sudden stall at an angle of attack of 10 degrees. These changes 
cause the peak value of power coefficient to be reduced slightly* 
and to occur at a hl^er tip speed ratio. Also* the no load point 
drops to about 8.6. For the catenary configuration the angle of 
attack variation along the blade is shown* on figure 7 for R ■ 4 
and R a 5* and at one angular location 9 » 0. As can be seen Hie 
angle of attack is almost constant over the majority of the blade and 
doesn't Increase appreciable until the last IS percent of the blade 
Is reached. However* the contribution of this part of the blade to 
the total power production is small. Shown on figure 8 is the 
variation of peak angle of attack with R at the maximum radius. 

From this it can be seen that at values of R ^ the peak angle of 
attack exceeds 10° consequently stall begins to set in. Since the 
angle of attack is essentially constant over the major part of the 
blade the whole blade stalls at about the same time thus causing 
the rapid reduction in the power coefficient curve as shown in figure 6. 

If a system could be built with no losses other than aerodynamic 
its starting characteristics would consist of an initial rotation 
caused by the differential drag on the two blades. However* the angles 
of attack experienced would be higher than the stall angle and limiting 
spin rate would be reached where the moment generated by the lift force 
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would juat balaaca that genera tad by the dreg force. To get the wlnd- 
nlll over thla point aona energy miat be added to the ayaten. Once 
the apln rate attaina a value high enough such that the net torque 
due to lift ia greater than that due to drag the aystera will spin up 
to the velocity ratio corresponding to no load condition. 

The computed results for the case n » 2 are also snown plotted on 
figure 9 versus the Induced velocity at the blades V. Included on this 
figure is a plot of the theoretical maximum curve obtained by setting 
the drag coefficient of the airfoil to zero. 

e 

As indicated previously no power is available until R reaches a 
value of about 2.5. As the spin rate Increases the power curve 
follows the trend of the theoretical maximum curve although it is 
displaced somewhat. As the spin rate increases to values greater than 
about 6 the drag losses becoote more Important since the angled of 
attack are decreasing and the contribution of the lift forces to the 
turning moment become of the sai^ order as the retarding effects of 
the drag forces. Consequently the actual curve departs from the 
theoretical maxlmuia by a large amount. 



23 


LINEAR THEORY FOR STSAIGRT BLADE COMFIOUBATIONS 


By restricting the possible windmill configurations to. 
those with straight blades and by assuming that the blade 
aerod 3 raamlc coefficients vary linearly with angle of attack, 
closed form analytical expressions for windmill performance 
parameters are obtained. These exptesslons are used to study 
the effects on windmill performance of variations in aspect ratio 
and solidity. 


Reduction to Analytical Form 

Consider the special case of a strcd^t blade configuration where 
r » 1.0 and 0 ■ 90°, both constant with respect to y. From figure 3 
we have 

V sin a = V cos 0 (33a) 

r 

and 


cos a “ R - V sin 0 

Also, if only small angles of attack are considered 

sin a 



VH.th these approximations the trust coefficient for one blade becomes 


(33b) 


(34a) 

(34b) 


C.p(0) 


% {[K, sin a cos a cos 0 + K, sin^ a sin 0] 
^ I X 



[sin a cos 0 - cos a sin 0]} 


V 

r 


2 


(35) 
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Substituting for yields 

CjO) - I { [ V (R - V sin Q) cos^ 0 + 

+ Cjj [V - R sin 0]} 
o 

and integrating over 0 an average thrust coefficient 
is obtained, 

K. TO + S| / [(R - V sin 0)2 

* 0 

s 

+ (V cos 0)2]^ (V - R sin 0) d ( 

Sl^llfying by assuming (R - V sin 0)» V cos 0 this 
nc K TO S 

^ / (R - V sin 0) (' 
^ 0 

or 

Oj. » |£RK V+3f£ VRCjj 
a o 

From momentum theory we have 

- 4 V (1 - V) 

^a 

therefore, equating these yields an expression for V 
velocity at the blade 

-if (Ki + 3 C„ ) 


y2 cos2 0 sin 0] 

(36) 

for n blades 
a 

) (37) 

becomes 

r - R sin 0) d 0 

(38) 


, the nondimensional 


V - 1 - 


(39) 
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Using eq. (9) vith r > 1 and 6 • 90°, th« aosent coefficient <^(6) 
can be written as 

2 Y 

CjjO) ” ^ 2 ; — c [ (Cj^ sin (* - Cjj cos a) 

^8 

+ c(.5 - d/c) (Cj^ cos a + Cjj sin a)l (^) 

but for a straight blate configuration A » 2 h 

8 

Using the relationships between « V in equation (33) 

yiel& * 

^ ^ —9 9 — _ _ — 

Cjj(0) “ 2 f ^ *^0® 0 + c (.5 - d/c) V cos 0 (R - V ain 0)] 


+ ~ (R - V sin 0)^ + cos^ 0 [c (.5 - d/c) V cos 0 

o 

- Cjj (R - V ain 9) ] 
o 


(41) 


Restricting this analysis to cases in which R - V sin 0 » V cos 0 
equation 41 can be integrated over 9 to yield the average moment 




nc Cp 

- ° (2R^ + V^) 


( 42 ) 


The power coefficient Cp is given by Cp = R , 

a 

substituting yields 
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C, - « ^ (1 - ^ 04 + 3 Cj ))* 

o 

s 

- (2R* + (1 - ^ + 3 ))^I 

o 

where 

1 - a^/v AR 

^1 “ *o + a /w AR^ 
o 


How for a atralglit blade 

AR - - - h/c 
c 


(43) 


(44) 


therefore. 


K 


1 


1 - a c/k h 
o 

1 - a c/"K h 
o 


(45) 


EquatioD (43) can be used Co evaluate the effects of nc, h, 

and on perfonance of wlndnllls . For the case where » 0 
o o 

the vd.ndi3llls performance reaches the theoretical maxliBum. This 

case along with a case using a representative value of are 

shown in figure 10. Both and R are shown as a ftmction of 

velocity ratio V. Also shown on this figure are the results 

of the numerical analysis using nonlinear aerodynamics. As would 

be expected the linear and nonlinear results are in fairly good 

agreement at the higher spin rates; however, comparing these with 

the limiting value given by the = 

*^o 


0 curve Indicates that 
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8lgnlflc«Dt i^provoKats in perforanoe could be realised if the 
airfoil drag could be reduced. For low spin rates the analytical 
and niaerical results are in poor agreenent due to the stall effects. 

The use of linear aerodynaatcs in the analytical solutira cai^es 

the drag tern to appear erroneously ssiall so that the analytical 

results approadi the theoretical narLaun for low spin rates. This 

is due to tiie fact that the airfoil experiences high angles of 

attadc and the contributicm of the lift forc» is orders of 

■agnitude greater than die retarding effect of the zero lift drag force. 

a 

Paraaeter Study 

The analytical soluticns can be used to evaluate the 
effects of solidity, and windmill aspect ratio on power coefficient. 

For snail values of solidity, which is proportional to nc, and for 
hi^ tip speed ratios these results are adequate to Identify trends; 
however, at low tip speed ratios the effects of nonlinear airfoil 
diaracteristlcs become dominant and the analytical results are 
inaccurate. Figure 11 shows the analytical results as a function 
of solidify at a fixed aspect ratio. Power coefficient is plotted 
versus tip speed ratio and as indicated on the figure peak 
>ifficien <7 Increases as solidity increases; however, the peak 
value occurs at a lower tip speed ratio and the range of tip 
speed ratios at which power can be supplied by the windmill varies 
inversely with solidity. Figure 12 shows the variation of power 
coefficient with it as tip speed ratio varies. As can be seen h 
does not have a major effect on windmill performance for the range 
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of values cited. In general aa h increases the aspect ratio of 
eadi blade Increases consequently its perforaance inproves. However, 
diK to structural considerations it is desirable to design a 
wln&lll with as snail a value of ii as possible witiiout drastically and 
inpalring perforaance. 

Hben a noKrlcal analysis based on the general theory is used 
to deteradne the effect of solidity on perfomance a major difference 
occurs for mall values of tip speed ratio. These results are 
shomt on figure 13. Since in this analysis the occurrence of 
airfoil stall is not d^endent on solidity the minimm speed ratio 
at which power is produced is essentially the sane for all 
conf Iguratlcms . Using the aerod 3 rnamic data for the NACA 0012 
airfoil this value of tip speed ratio is about 2.5. The efficiency 
increases rapidly as tip speed ratio increases from this minimum 
value. Hwever, for seme values of solidity the peak efficiency 
would occur at tip speed ratios less than this minimum value were 
it not for stall effects. Consequently, with stall effects 
included, these peak values are never achieved. If a windmill could 
be designed using a blade which did not stall or stalled at a 
hi^?r angle of attack it might be possible to achieved the hl^er 
performance indicated by the analytical solution. However, it 
remains questionable if one would want a design which responded 
so radically to changes in wind velocity. 
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UOHFIGORAXION EFFECTS Oil AERWYNAMIC EFFICIEHCT 

Considering only aerodynaaic efficiency the optimum VAW design 
as depicted in figure 14a would consist of straight blades attached 
to the center shaft by low drag struts. The wlndaill aspect ratio 
■ust be large enough to adniaize losses due to blade aspect ratio. 
Unfortunately such a design leads to severe structural llaitations 
due to the hi^ bending loads which develop in the blade. It is for 
this reason that curved blades have been selected for designs currently 
In operation. With a curved blade the*bending loads are greatly 
reduced and the centrifugal loads are reacted as almost pure tension 
loads in the blades. However, it is possible to desi^ a straight 
bladed VAW which has the same structural cap^ility as a curved blade 
VAW. This can be accomplished by supporting the blade at appropriate 
intervals as shown in figure 14c. Secondly, by placing the support 
struts at an angle other than 90 degrees from the axis of rotation 
they will produce enough lift force such that their lift induced torque 
is greater than the drag induced torque ar.d they in effect contribute 
to the performance of the windmill. 

The results of such a design are shown in figure 15 along with a 
comparison of the catenary VAW performance. For this configuration 
the support struts and the primary blade are an NACA 0012 airfoil. The 
strut to blade chord ratio is 1/2. As the data Indicates the peak 
power Coefficient for the chevron design is about .43 compared to the 
value of .37 for the catenary. However, the catenary can produce 
power over a wider range of spin ratios than the chevron. 
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Ho has be«n Mdo to optlal&o tbe dievroo dMlgo and one 

would expect that for higjber valuea of h/r soae ieproveaent in 
perfotaanoe could be expected. Also, the hi|^er tiie angle of the 
support struts the greater will be their contribution. It is also 
possible to select a support strut mgle sudi that the angle of 
attack at soae point on the strut is the sane as the angle of attadc 
at the blade. 




CCmCLUSIONS 


The results from the analysis presented herein Indicate that a 
nanerlcal solution using non-linear aerodynamic data and assuming a 
uniform Induced velocity over the coiqilete volume swept by the windmill 
can be used to predict performance of vertical axis windmills. The 
major parameters affecting windmill performance are the blade shape, 
blade aspect ratio and downwash effects. The effects of curved flow 
are insignificant providing the ratio of blade chord to diameter is saiall.. 
The use of linear aerodynamics allows t)ie defining torque equatlcm to be 
integrated; however, these analytic solutions are Incorrect at lower values 
of spin parameter where the aerodynamics are nonlinear. They are useful 
for assessing trends at higher spin parameter values. The performance of a 
VAW can be improved by using straight rather than curved blades provided 
the blade is supported in the proper manner. However, optimization of 
any VAH design will ultimately be based on the costs per Installed unit of 
power consequently not only aerodynamic but also structural considerations 
must be taken into account. 
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APPENDIX 


WINDIfILL PERPORHANCE PREDICTION COMPUTER PROGRAM 


Order of Input 

1. Header Card, Colunos 1 through 80. 


2. Maaellst $NAM1 

including: 


Fortran 

Sviri>ol 

Engineering 

Synbol 

Definition 

AR 

AR 

Aspect ratio of blade 

ARS 

AR 

s 

Aspect ratio of lifting struts 

AZ 

a 

z 

Section lift curve slope 

BETA(I) 

6 

Array of local slope of blade 
for each YBAR(I) 

BETS 


Angle of lifting struts 

cbar 

c 

Nondin^nsional blade chord 

CBARS 

c 

s 

Chord of lifting and nonlift lug 
struts 

CD(I) 

s 

Array of Section drag coefficient 
of blade and struts for each 6AH(I) 

CDl 


Zero lift drag oefficient of 
nonlifting struts 

CKl 

CKl 

Number of nonlifting struts 

CL(I) 

s 

Array of Section lift coefficient 
of blade and struts for each GAM(I) 

CMC4(I) 

X/4 

Quarter chord pitching mon^nt 
coefficient 

DY 

DY 

IncrezDent used co calculate y 
coordinate array for blade, from 
YMIN to YMAX* 

DYS 

DY 

8 

Increment used to calculate Y 
coordinate array for struts from 
YSMIN to YSMAX. 

GAM(I) 

Y 

Array of angle of attack, 
where I » 1, NGAM. 
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Fortran 

Engineering 


Symbol , 

Syidjol 

Definition 

ITHPENT 


Total number of THPRNT values* 
maximum is 10. 

N 

n 

Number of blades 

NGAM 


Total nund>er of GAM values* 
maximum is 100. 

NS 

n 

8 

Number of lifting struts* 
if none set = 0. 

NTHETA 


Total number in THETA 
array* maximum Is 200. 

NVR 


Total number in VRRY 
array* maximum is 10. 

NY 


Total number in YEAR array* 
maximum is 50. 

RBAR(I) 

r 

Array of windmill blade raditis 
for each YBAR(I). 

TC 

T 

Blade airfoil thickness ratio. 

THETA(l) 

e 

Array of windmill angular 
coordinate* degrees* where I ■ 1* 
NTHETA. 

THPENId) 


Array of THETA values at whidi 

the lift* drag- and moment values 

can be printed, where 1=1* ITHPRNT. 

VRRY(I) 

R 

Array of specified velocity ratios 
at which the windmill is to be 
analyzed* where 1=1, NVR. 

YBAR(I) 

y 

Array of coordinates along the 
axis of rotation of the windmill* 
where 1=1, NY. 

YMAX 

Y 

max 

Maximum Y coordinate of blade. 

YMIN 

Y . 
min 

Minimum Y coordinate of blade. 

YSMAX 

Y 

smax 

Maximum y coordinate of lifting struts 

YSMIN 

Y . 
smln 

Minimum y coordinate of lifting struts 
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Output Parameters 


Output 

Svnbol 

Engineering 

SyiN>ol 

Definition 

BARCI 

S: 

a 

Windmill thrust coefficient averaged 
over a revolution 

CD 


Airfoil section drag coefficient 

CL 


Airfoil section lift coefficient 

CM(TBEIA) 


Bla<te Boment coefficient for angular 
•location G 

OlA 

Si 

a 

Average windsdll moment coefficient 

CMi^ 

Si 

as 

Average moment coefficient due to struts 

CMC4 

^“c/A 

Airfoil section quarter chord mcmient 
coefficient 

QIS (THETA) 

s 

Strut moment coefficient for 
angular location 0 

cm 

Cm 

Total moment coefficient for n blades 
at angular location G 

CUTS 


Total moment coefficient for all 
struts at angular location G 

CPA 

a 

Pever coefficient due to lifting 
blades and struts 

CPL 


Power coefficient loss due to 
nonlifting struts 

CPN 

s 

n 

Net power coefficient 

D/C 

d/c 

Location of center of pressure 


measured chordwlse from the airfoil 
leading edge in fractions of the 
chord length 
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Output 

Svid)ol 

Engineering 

Svabol 

Definition 

DEL CM 


Incremental moment coefficient due 
to a blade e].ement 

GAMMA 

Y 

Local blade angle of attack 

OMEGA 

0> 

Windmill rotational speed 

EMAX 

^max 

Maximum windmill radius 

SBAR 

S 

Nondlmenslonal windmill projected area 
in a vertical plane 

THETA 

0 

Angle of windmill rotation measured 
from the free wind ulrection 

VINF 


Wind stream speed at a large distance 
from the windmill 

VR 

V 

Wind stream speed at the windmill 
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PROGRAM FLOW DIAGRAM 
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2 



CALCULATE FIRST ESTIMATE 
OF VR, VELOCITY AT WINDMILL 
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Table I. Aerodynamic Characteristic of the 

NACA 0012 Airfoil at Reynolds Numbers 
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Figure 1. The Vertical Axis Windmill. 






Windmill blade viewed from above. 



View B>B 


Figure 3. - Velocity and force vector components. 




Local Coordinates for Flow Curvature 
Analysis. 






Figure 5. - Side View of the Upper Half of a Blade. 







- Variation of Angle of Attack Along 
the Blade - Catenary VAW. 


FIGURE 7 
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FIGURE 8 - Maximum Angle of Attack Envelope for 
the Catenary VAW as a Function of R. 
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FIGURE 9 - Performance Characteristics of the 
Catena ry Wi ndmi 1 1 . 



Analytic Solution Numerical Solution 
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FIGURE 10 - Performance Characteristics of a 
Straight Blade VAW. 
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STRUTS 




b) TYPICAL LOADING FOR STRAIGHT AND CURVED BLADE VAW. 



c) CHEVRON STRAIGHT BLADE VAW CONCEPT. 


FIGURE 14 


Straight and Curved Blade VAW Designs. 







